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ABSTRACT 

Gamma-ray burst (GRB) host galaxies have been studied extensively in optical photometry and spectroscopy. 
Here we present the first mid-infrared spectrum of a GRB host, HG 031203. It is one of the nearest GRB hosts 
at z = 0. 1055, allowing both low and high-resolution spectroscopy with Spitzer-IRS. Medium resolution UV-to- 
Zf-band spectroscopy with the X-shooter spectrograph on the VLT is also presented, along with Spitzer IRAC 
and MIPS photometry, as well as radio and sub-mm observations. These data allow us to construct a UV-to- 
radio spectral energy distribution with almost complete spectroscopic coverage from 0.3-35 /im of a GRB host 
galaxy for the first time, potentially valuable as a template for future model comparisons. The IRS spectra 
show strong, high-ionisation fine structure line emission indicative of a hard radiation field in the galaxy - in 
particular the [S IV]/[S III] and [Ne III]/[Ne II] ratios - suggestive of strong ongoing star-formation and a very 
young stellar population. The absence of any PAH emission supports these conclusions, as does the probable 
hot peak dust temperature, making HG 03 1203 similar to the prototypical blue compact dwarf (BCD) galaxy, 
IIZw40. The selection of HG 03 1203 via the presence of a GRB suggests that it might be a useful analogue 
of very young star-forming galaxies in the early universe, and hints that local BCDs may be used as more 
reliable analogues of star-formation in the early universe than typical local starbursts. We look at the current 
debate on the ages of the dominant stellar populations in z ~ 7 and z ~ 8 galaxies in this context. The nebular 
line emission is so strong in HG 03 1203, that at z ~ 7, it can reproduce the spectral energy distributions of 
z-band dropout galaxies with elevated IRAC 3.6 and 4.5 /im fluxes without the need to invoke a 4000 A break. 
Indeed, photometry of HG 03 1203 shows elevation of the broadband V -magnitude at a level similar to the IRAC 
elevation in stacked z-band dropouts, solely due to its strong [O III] line emission. 

Subject headings: gamma-ray burst: general - gamma-ray burst: individual: GRB03 1203 - early universe - 
dark ages, reionisation, first stars - galaxies: dwarf - galaxies: ISM 



1. INTRODUCTION 



Long-duration gamma-ray bursts (GRBs) trace the f orma 
tion of massive sta rs ( Galama et al.| 199 8 1 Hjorth et al .|2003 
|Stanek et al.|2003| |Bloom et al.|2002| i~They occur across a 



vast r ange of cosmic distances from z = 0.0085 (Galama et al. 
1998] > to z = 8.2 dTanvir etal.|2009ll, an d with a very high mean 
redshiift ((z) = 2.2 Fynbo et al. 2009). Their hosts therefore 
provide a unique way to select star-forming galaxies because, 
unlike other selection techniques, their selection is unrelated 
to the emission properties of the galaxies themselves. 
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Studies of GRB hosts have so far suggested that they 
are star-forming galaxies with relatively high s pecific star- 
formation rates and young stellar populations (IC hristensen 
etlLl[2004l ICas tro Cero n et al 1|2006| |Savaglio et aLpoTF^ 
Svensson et al. 2010; Chary et al. 2002). However, it is worth 
bearing in mind that such general statements are still difficult 
to make in an absolute sense because of the potential bias in- 
troduced in sample selection because most GRB hosts have 
been selected up to now on the basis of localisations depen- 
dent on optical afterglows. Some attempts to address this es- 
sential question are b eing made with X-ray-selected samples 
( |Malesani et al.|2009] >. 

Optical spectroscopy in absorption using the afterglow (e.g. 
Fynbo et al .|2009 1, or occasionally in emissi on from the host 
galaxy itself w here it i s bright enough (e.g. Wiersema et al. 
|2007| Levesqu e et al.||2010| l have indicated metallicities in 
GRB host galaxies that are relatively low, but not atypical 
of the population forming high-ma ss sta rs, at least at z > 2 
( |Fynbo et al.l|2006l [Caiura et al.|[2009> . One of the best- 
studied cases so far, because of its pr oximity at z = 0.1055, i s 
the host of GR B 031203 (HG031203,IProchaska et al"] 2004). 
These studies (jP rochask a et alj|2004j |Margutti et al'j|2007l 
Levesque et al. 2010||Han et al.|2010| i have shown HG 03 1203 
to be a relatively low metallicity system with a fairly high 
star-formation rate and a young stellar population. As with all 
optical studies, such analyses are dominated by the low-dust 
regions and may give an incomplete impression of the galaxy 
since a significant fraction of the stellar emission, particularly 
from young stars, may be re-processed by dust and emerge in 
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the mid- and far-infrared (IR). 

To date, all spectroscopic studies of GRB host galaxies have 
been at optical or near-IR wavelengths. Here we present the 
first spectroscopic study of a GRB host galaxy in the mid-IR 
(5-40 /im) using the Spitzer Space Telescope (Wern er et al. 
12055) Infrared Spectrograph (IRS |Houck et al.||2004) . We 
also present an almost continuous spectral energy distribution 
(SED) for this galaxy from the ultraviolet through to the mid- 
IR (0.35^-0 /im), with detections and upper limits at far-IR, 
sub-mm and radio wavelengths, the first such SED for a GRB 
host. 

Uncertainties quoted are at the 68% confidence level for 
one interesting parameter unless otherwise stated. A cosmol- 
ogy where Hq = 72 km s" 1 Mpc" 1 , £7 a = 0.73 and il m = 0.27 is 
assumed throughout. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. X-shooter 

HG 03 1203 was observed during the commissio ning of the 
second gen eration VLT instrument, X-shooter ( |D'Odorico| 
|et al.|[ 2006) on 17 March 2009. X-shooter is a single object 
echelle spectrograph consisting of three arms which cover si- 
multaneously the spectral range 300-2480 nm. The slit widths 
chosen during the observations were 1"0 in the UVB arm 
(AA = 300-550 nm), and 0"9 in the VIS (AA = 550-1015 nm) 
and NIR arms (AA = 1000-2480 nm). This instrument setup 
gives resolutions in the three arms of R = 5100, 8800, and 
5600, respectively. The slit width matched the seeing during 
the observation. During the observations the slit was aligned 
along the parallactic angle. We obtained four 20 minute in- 
tegrations on target, with offsets along the slit in a classic 
ABBA observing pattern. The UVB detector was binned by a 
factor of two in the spectral direction when the detectors were 
read out. With a 1" slit width in the UVB arm, the sampling 
for the full-width half-maximum of a line is 6 pixels, so a bet- 
ter signal-to-noise ratio per pixel is obtained in binned data. 

Data reductio n was performed with a preliminary version of 
the pipeline (Goldon i et al.|2006) . Tracing the individual or- 
ders, wavelength calibration, and flat fielding were done using 
calibration frames obtained during the commissioning run. In 
the UVB and VIS data, subtraction of the sky background 
used a routine described in Kelson] d2003), while the back- 
ground in the NIR data was subtracted using the adjacent sci- 
ence exposures. Finally the individual orders were extracted 
and merged using a weighted mean combination scheme, re- 
sulting in a two-dimensional spectrum for each arm. Similar 
data reduction procedures were done for the HST white dwarf 
spectrophotometric standard star GD71 which was used to 
flux calibrate the science data, and for an observation of 
the 08V star Hipparcos 69892, which provided a reference 
smooth spectrum for dividing out telluric absorption lines 
from the earth's atmosphere. 

Further analysis of the data was performed with conven- 
tional IRAF routines. The four two-dimensional spectra were 
combined, and one-dimensional spectra were extracted using 
an aperture which matched the three arms. In the flux cali- 
bration of the observation, we included a correction for the 
atmospheric extinction using an extinction curve appropriate 
for Paranal. 

2.2. Spitzer 

The Spitzer photometric data on HG 031203 were obtained 
with IRAC ( [Fazio et al.|2005) at 3.6 and 5.8 /zm with a 30 s 



Table 1 

Spitzer-IRS observations of the host galaxy of GRB 031203 



Module 


SH 


LH 


SLl 


SL2 


LLl 


LL2 


Integration Time (s) 


3656.9 


1828.5 


609.5 


1828.5 


314.6 


943.7 



frame time, and a total integration time of 360 s and MIPS 
(Riek e~et al.|2004) at 24 /im also with 30 s integration times 
and a total observation time of ~ 1350 s. The data, under pro- 
gram 2 0370, were reduced as outlined in Castro Ceron et al. 
(|20T0j. 

The target was also observed spectroscopically, in star- 
ing mode using all four IRS modules; short-low (SL, 5.2- 

14.5 /im), long-low (LL, 14.0-38.0 /im), short-high (SH, 9.9- 

19.6 /im), and long-high (LH, 18.7-37.2 /im). The resolution 
of the low-resolution modules is ^60-127, while that of the 
high-resolution modules is ^600. Observations were carried 
out on 28 May 2005 (see Table[jl 

Data from the short-low and snort-high modules were pre- 
processed with the Spitzer Science Centre (SSC) data reduc- 
tion pipeline version 15.3, while data from the long-low and 
long-high modules were processed with version 17.2. 

For the high-resolution modules, spectra were extracted 
from each individual basic calibrated data (BCD) file using 
the Spitzer IRS Custom Extraction (SPICE) software. Prior to 
extraction of the spectra, data files were individually cleaned 
of rogue pixels using the IDL program IRSCLEAN version 
1.9. The rogue pixel masks used for cleaning were a com- 
bination of the default campaign masks provided by the SSC 
and masks generated automatically by running IRSCLEAN's 
rogue pixel identification algorithm. 

Individual spectra were subsequently combined with a 
clipped mean algorithm using the Spectral Modelling, Anal- 
ysis, and Reduction Tool (SMART) developed by the IRS 
team (Higdon et al. 2004). The edges of each order were 
trimmed based on the wavelength calibration ranges provided 
in Table 5.1 of the IRS Data Handbook version 3.1. Since 
no sky measurements were taken, the contribution of the sky 
emission was not subtracted from the high-resolution spectra. 
No subtracting the sky background has a significant effect on 
the level of the continuum measured with the high-resolution 
modules. However, we use the high-resolution data only to 
extract emission line fluxes, not for the SED work, so that 
a continuum offset is not a significant issue. The fluxes ob- 
tained with the high-resolution modules are consistent with 
those retrieved from the low resolution spectra. 

For the low-resolution data, we made use of pipeline co- 
added, sky-subtracted products. These files result from co- 
adding the data from each nod position, then subtracting each 
nod position from the other in order to subtract the back- 
ground. The spectra were cleaned for rogue pixels and op- 
timally extracted to produce the final 1-D spectra. 

The source was also observed with MIPS in SED mode. 
The data were re duced in a standard way, as described in | van 
|Loon et al.| | |20l0) . 



2.3. Sub-mm-radio 

We obtained radio continuum observations of HG 03 1203 
on 26 January 2008 using the Australia Telescope Com- 
pact Array (ATCA) in configuration 6A, with antennae po- 
sitioned on both east-west and north-south tracks, and base- 
lines of 60^-500 m. Simultaneous observations were made 
at 20 cm (1.39 GHz) and 13 cm (2.37 GHz), with a bandwidth 
of 128 MHz at each frequency. A total of 7 hr of on-source 
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Rest Wavelength (\im) 
11 12 13 15 




Observed Wavelength (jj,m) 

Figure 1. Low-resolution IRS spectrum of the host galaxy of GRB 031203. Orders 1 and 2 of the SL (left) and LL (right) modules. Thin red lines indicate 
uncertainties on the spectrum. Strong fine structure lines typical of a young star-forming galaxy are detected (see Fig. l2]for close-ups of the lines as observed 
with the high-resolution modules). Very high-ionisation fine structure lines, indicative of AG N activity (e.g. [ Ne V] 14.3 /^m or [O IV] 25.89 (im, |Tommasin| 
|et al.|20T0) , are not found. The IRS spectrum of the blue compact dwarf galaxy IIZw40 from |Wu et al.|(2~006fr is plotted (dashed line), scaled and oftset, tor 
comparison purposes. 



data were obtained. Calibrator sources PKSB1934-638 and 
PKS B0826-373 were utilized to set the absolute flux cali- 
bration of the array and to calibrate phases and gains, re- 
spectively. Data reduction and analysis was done using the 
MIRIAD package ( |Sault & Killeen|2004| l. The final synthe- 
sized beam sizes for 20 and 13 cm images were 8 "5 x 3. "4 and 
6"3 x 2"3, respectively, with root-mean-square (rms) values 
of 46 and 37 fi]y beam" 1 . The flux density of HG 03 1203 was 
estimated by fitting a two-dimensional Gaussian function to 
the data with its centroid, size, and orientation as free param- 
eters. 

We obtained sub millimeter (870 /im) observations on 13- 
18 August 2008 using the Large Apex BOlometer CAmera 
(LABOCA; |Siringo et~aLl|2"00 9) mou nted at the Atacama 
Pathfinder Experiment (APEX; Gi isten et al.||2006| l. A total 
of 7 hr of on-source data were obtained . The we ather condi- 
tion varied between 0.3-1 .4 mm of precipitable water vapour. 
Data red uction and analys is was done using the miniCRUSH 



Table 2 

Spitzer, APEX, and ATCA photometry of HG 03 1203 



Instr. 


IRAC 


MIPS 




APEX 


ATCA 






(mJy) 




(mJy) 




Band 


3.6 /xm 5.8 /im 


24 /im 


SED 


870 fim 


13 cm 21cm 


Flux 


193 ±2 356 ±13 


11.3±0.4 


<40 


< 12 


191 ±37 254 ±46 



Table 3 

Principal emission lines in the mid-IR spectrum of HG 03 1203 



package ( |Kovacs|2 008 [ 



We used the 'deep' option that re- 
sults in the best signal-to-noise ratio for faint, point-like ob- 
jects. The beam size for the final image was 19 "5 and the rms 
was 2.2 mJy beam" 1 . The target was not detected. 

3. RESULTS 

Spitzer, APEX and ATCA each provide photometry or pho- 
tometric limits at long wavelengths, allowing us to determine 
the extinction-free emission from HG 031203. These data are 
presented in Table [2] The Spitzer photometric detections are 
consistent with the flux-calibrated IRS spectra at the same 
wavelengths. 

3.1. Spitzer spectroscopy 



Line ID 


Wavelength (pm) 


Flux (10 15 erg cms ') 


[S IV] 


10.51 


8.5 ± 1.5 


[Ne II] 


12.81 


0.7 ±0.2 


[Ne III] 


15.56 


10.6± 1.1 


[S III] 


18.71 


5.2 ±0.7 



The most striking features of the IRS spectra of HG 03 1203 

lines, [S IV] 10.51 /xm, 
and [Ne II] 12.81 /im 



(Fig. [Til are clearly the strong forbidden lines, [S IV] 10.51 /im, 
[S IIIT18.71 /Lim, [Ne m] 15.55 ^m, 
(Fig.[|andTable[3]>. 

Comparing our spectrum to spectra of other star-forming 
galaxies ( |Wu et al.|2006[l, it is immediately apparent that the 
broad features at 6.2, 7.7, 8.6 and 11.2, 12.7, and 16.4 /im, 
believed to be a ssociated with polycyclic aromatic hydrocar- 
bons (PAHs, |Tielens||2008| l, are, somewhat surprisingly, not 
detected. The lack of PAH emission could be indicative of a 
very strong continuum from an active galactic nucleus (AGN) 
diluting the PAH features, a discriminant noticed with ISO 



http : //www . submm. caltech . edu/~ share /crush/ 



(e.g. |Lutz et al.|[T998 Laureijs et al.|2000 l. The presence of 
an AGN is claimed by JLevesque et al. (2010), however, even 
at first glance this seems unlikely since the PAH emission ap- 
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Figure 2. The locations of the most prominent mid-IR fine structure lines 
typically found in star-forming galaxies from the high-resolution, SH and 
LH, modules. The grey lines indicate the associated uncertainty at each wave- 
length. The broad feature at 10.2 fj,m is undetected in the SL spectrum and is 
related to an unremoved artefact in our SH data. 

pears to be simply absent, rather than diluted by a powerful 
continuum. We address the other lines of evidence that essen- 
tially rule out a significant AGN contribution in this galaxy 
below. The other possibility is that the galaxy is a power- 
ful starburst. In such galaxies, it is well-known that the PAH 
flux is anti-correlated with the hardness of the radiation field, 
either because young star-forming populations do not form 
PAHs, or b ecause the hard radiation field destroys them ( |Wu| 
|et al.|20 06). This seems consistent with our findings from the 
mid-IR and optical-near-IR spectra below. 

The ratios of the forbidden lines, associated with the ex- 
citation state of the gas in the interstellar medium of the 
galaxy, indicate the hardness of the radiation field. In Fig. [3] 
we show the [S IV]/[S III] flux ratio (1.6tJ|) as a function 
of [Ne III]/[Ne II] (16.3 ±4.5), and compare HG 031203 to 
samples of starburst galaxies and blue compact dwarf galax- 
ies (BCDs). The spectrum shows no indication of absorp- 
tion at 9.1 pm. from silicates, neither does it show any hint 
of the higher excitation lines that might be expect ed from an 
AGN such as [Ne V] 14.3 pm or [O IV] 25.9/zm ( |Tommasin| 
[eTlLl[20T0| . All of these facts point to HG 031203 being a 
purely star-forming galaxy, but with an exceptionally hard ra- 
diation field, in dicating a very y oung stellar population, prob- 
ably < lOMyr ( |Madden|2000| ). The fact that these ratios are 
observed in the MIR suggests that the galaxy has essentially 
the same character throughout, i.e. that it is dominated in the 
MIR by the same young stellar population that dominates in 
the optical (see below). 

3.2. X-shooter spectroscopy 

The X-shooter spectra are very information-rich (Fig. [4}, 
containing many strong emission lines with a resolving power 
high enough to look in some detail at the interstellar medium 
(ISM) in this galaxy. However in this paper we will address 
ourselves primarily to the continuum, with regard to the anal- 
ysis of the SED of the galaxy in order to determine the domi- 
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Figur e 3. Line ratios of star-fo rming galaxies. The starburst galaxy sam- 
ple of |Bernard-Salas et al. 1 2009 1 are plotted as circles. A sample of blue 
compact dwart galaxies, which have younger dominant stellar populations, 
higher specific star-formation rates, and harder radiation fields even than the 
starburst galaxies, is also plotted (open squares, Wu et al. 2006 1. HG 03 1203 
(filled square) is among the most extreme of the blue compact dwarf galaxies, 
indicating that its dominant stellar population is extremely young. The fact 
that the diagnostic lines are in the mid-IR ensures that the diagnostics are not 
significantly affected by dust obscuration. 



nant emission processes and the overall nature of the ISM. In 
this regard, we have also determined the ratios of the fluxes of 
nebular lines, to determine the nature of the source powering 
the optical emission: [O II], H(3, [O III], [S II], [N II], and [O I]. 
The relevant line ratios are: log([0 III]/H/3) = 0.77 ± 0.02, 
log([S E]/Ha) = -1 .23 ± 0.02, log([N II]/Ha) = -1 .35 ± 0.02, 
log([0 q/Ha) = -2.02^$, and log([0 III]/[0 II]) = 0.78;$$ . 

These high-quality ratios resolve the debate about the pres- 
ence of an AGN in HG 03 1203. While the line ratios do in- 
dicate a very hard radiation field in the ISM of this galaxy, 
making it among the most extreme starburst galaxies con- 
sistent with the MIR line ratios, they unequivocally locate it 
in the star-forming region of the Baldwin-Phillips-Terlevich 



(BPT, Bald win et al.||198T) excitation diagrams (see Kewley 
|et al.|2 006 and referen ces therein), contrary to a recent claim 
( |Levesque et al.|2010|l, and consistent with previous findings 
of a pure sta rburst (P rochaska et al.|2004||Margutti et al.|2007| 
|Han et al.|2010[ ). Our similar finding from the mid-IR spec- 
trum above on the absence of AGN characteristics, indicates 
that there is no very different obscured component. 

The presence of a break in the spectrum at 4 000 A is an 
indication of the age of the stellar population ( Bru zual "&] 
|Charlot|[2003] l. We quantified the break in HG 03 1203 us- 
ing the definition of D4000 of |Balogh et al.| ( |1999") , finding 
D4000 = 0.93 ±0.05. This effectively means thatThere is no 
4000 A break, and indicates a continuum dominated by a very 
young stellar population, probably younger than ~ lOMyr, 



with a high specific star- formation rate (Brinchmann et al. 
|2lX^|GaiTazzi et al.|2005| ). 

3.3. The spectral energy distribution ofHG 031203 
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Figure 4. The X-shooter spectrum of HG 03 1203. Strong nebular emission dominates the spectrum. 



The combination of 0.3-2.5 /im spectroscopy from X- 
shooter and 5-35 /im spectroscopy from the IRS, combined 
with limits in the FIR and sub-mm and detections at radio 
wavelengths, allow us to get a fairly complete picture of the 
properties of HG 03 1203 via its SED (Fig. [5}. 

3.3.1. Foreground dust extinction 

The first outstanding question is the Galactic dust obscu- 
ration toward this galaxy. The Balmer decrement, which, 
due to the low redshift, essentially gives us the to tal redden- 
ing along th e line of sight, yields E{B-V) =1.17 (Prochaska 
|et al. [2004): it is not entirely clear, however, what fraction 
of this extinction is in the Galaxy and what fraction is in 
HG 031203. Before discussing the foreground extinction, it 
should be noted that our main results for HG 03 1203 are in- 
sensitive to the precise value of the extinction fraction in our 
Galaxy. The mid-IR emission is largely unaffected by this 
level of extinction. Our primary age indicator for the galaxy 
population, D4000, is also insensitive to the precise extinction 
correction applied, as are the most important optical line ratios 
since they lie very cl ose in wavelen gth space (e.g. O III/H/3). 

The dust maps of ISchlegel eTaTI ( fl998] l yield E(B-V) = 
1 .04 in the Galaxy alone in this directi on and this value has 
been used by several groups f Thomsen et al. 200 4} |Gal-Yam| 
[eTaLl[2004l |Malesani et aLl|2004| |Cobb et al.||2004| i. While 
these measurements seem to imply that the correction fac- 
tor to obtain the intrin sic galaxy spectrum is E(B- V) = 1.04, 
Prochaska et al.|(|2004| and Margutti etaTT] <|2007)> in fact ar- 



along this sightline and suggest the foreground Galactic ex- 
tinction is E(B-V) = 0.7-0.8, with E(B-V) ~ 0.3- 0.4 resid- 
ing in HG 03 1 203 . Recently, ISchlafly et al] ( p010) l have in- 
dicated that the global normalisation for E(B- V) in Schlegel 
et al. (|1998|l may be overestimated by ~ 14%. There are, 



however, significant variations in the normalisation over the 
sky (~ 20%): in particular, in this direction in longitude, 
the normalisation is somewhat higher, but for high E(B-V) 
sightlines, the normalisation is somewhat lower. Furthermore, 
for very high extinction sightlines, the extinction-estimation 
method used is not effective. This leaves the situation some- 
what confused. 

Another line of attack to address the level of the Galactic 
extinction is to use the Leiden-Argentine-Bonn Galactic H I 
maps, which i ndicate a column den sity of 5.0 x 10 21 cm" 2 in 
this direction (Kalberla et al. 2005 ). The gas to dust corre- 
lation of Predehl & Schmitt ( 1995 ) and more recently |"Vuong| 



etal.(2003), suggest an E(B-V)df 0.91 and 1.03 respectively, 



assuming standard ISM abundances; but the measurement is 
abundance sensitive and may be lower by 20% for a differ- 
ent choice of metallicity ( |Vuong et al.|20 03 ). The situation is 
therefore still complex. 

However, in fitting the SED of HG 03 1203, we find that the 
total star-formation rate (SFR) is strongly inconsistent (larger 
than a factor of two) with the Ha-derived SFR if we assume 
a low Galactic extinction (E(B - V)=0.78). We req uire an ex- 
tinction close to the value of Schlegel et al. ( 1998) to recon- 
cile the SED-derived SFR with the SFR derived from Ha by 



gue that the value from |Schlegel et al.| ( |1998| ) is overestimated 



Prochaska et al.| ( |2004"] ) and Margu tti et"aT] (|2007). However, 
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the presence of a significant MIR dust emission in HG 03 1203 
suggests that the UV light suffers some dust extinction in the 
host, and our modelling of the full SED suggests Ay of a 
few tenths of magnitudes (see below). This is consistent with 
the observed Balmer decrement since it adds relatively little 
(E{B-V ) ~ 0. 1) to the overall reddening. On the basis of these 
arguments, we simply use E(B-V) = 1 .04, the value from the 



Schlegel et al. ( 1998) dust maps, to correct for Galactic fore- 
ground extinction in our analysis of the HG 03 1203 SED. 

3.3.2. Comparison to active blue compact dwarf galaxies 

The strong nebular emission lines, extreme hard line ratios, 
lack of PAH emission features, and probable hot dust tem- 
perature, suggest a close similarity between HG 03 1203 and 
blue compact dwarf ga laxies. This suggestion was made by 
|Sollerman et al.| ( f2005[ ) when analysing the optical properties 
of low redshift S5RjRB host galaxies, including HG 03 1203. 
In particular, HG 03 1203 is similar to those cl assified as being 
in an "active" mode of star-formation (e.g. |Hunt et al.|20Q5| >. 
The IR line ratios in fact, suggest that HG 03 1203 is among 
the more extreme BCDs in terms of the hardness of its radi- 
ation field (Fig. [3]). A striking similarity exists in the mid-IR 
emission of HG 03 1203 and II Zw 40, which ha s a very sim- 
ilar metallicity and line ratios ( |Wu et alT||2008) l and a nearly 
identic al IRS sp ectrum in terms of slope and strength of fea- 
tures (|Wu et al.|2006). We have plotted the SED of IIZw40 
froml Hunt et al.| ( |2005) on top of the SED of HG 03 1203 in 
Fig. [5 scaled by the ratio of the [Ne III] 15.6 /im line fluxes. It 
is clear that it is only in the mid-IR that the spectral similarity 
holds. The UV/optical continuum of II Zw 40 is considerably 
bluer; its IR peak is brighter than is possible for HG 03 1203, 
and is likely somewhat cooler. Both facts might be explained 
by a distribution of dust in II Zw 40 that is different, possibly 
distributed more uniformly, and less densely located in star- 
forming regions than is the case in HG 03 1203. 

3.3.3. Modelling the SED 

We attempted to fully reproduce the UV-radio SED o f 
HG 031203 using the GRASIL code dSilva et al.| [19981. 



The best-fitting SED model (see Fig. |5) was obtain ed wi 



th 



the fitting procedure outlined in MichaTows ki et aL] ( 2010a) 
using the templates of Iglesias-Paramo et al. ( |2007| l and 
Michalows ki et aL] ( |2008| l with the PAH emission removed 
and restricting the time of the total galaxy evolution to 
200 Myr. This attempt failed to reproduce the SED on two 
points. First the observed 3.6 flux was too high with re- 
spect to any of the models. This excess is well-known in 
dwarf galaxies and is speculated to be due to a hot dust com- 
ponent, strong Bra em ission, or nebular continuum emission 
(Smith & Hancock 2009). Given the effect of the very strong 
nebular emission lines on other photometric datapoints (e.g. 
the elevated V-band magnitude, Fig. Bl), line emission from 
Bra may be sufficient to explain the observed 3.6 /im excess 
here. The second point of failure relates to the 4000 A break, 
where none of the templates could reproduce the lack of a 
break in the spectrum HG 03 1203, in spite of the restriction to 
younger stellar templates. 

Finally, the dust mass is not very well-constrained with this 
dataset, since we do not quite cover the peak of the thermal 
emission component. An indication of a turnover in the IRS 
spectra at the longest wavelengths (Fig. [5) suggests the peak 
dust temperature is v ery high, as observed in dwarf galaxies 
(Galametz et al. 2009), perhaps hotter eve n than the 45 K sug- 
gested for sub-mm-detected GRB hosts (Michalowski et aL] 



2008 ). However our constraints at longer wavelengths do not 
allow us to exclude a large cool dust component. 

4. DISCUSSION 

We have found that HG 03 1203 is an extreme star-forming 
galaxy, similar to, but considerably more luminous than, 
nearby blue compact dwarf galaxies. A very hard radiation 
field is found from both optical and mid-infrared line ratios, 
and is consistent with a very young star-forming population. 
In spite of its blue optical continuum, approximately as much 
energy is reradiated by dust as comes directly from stars. 
This suggests a strongly non-uniform dust distribution in the 
galaxy. 

The host galaxy of GRB 031203 is among the brighter GRB 
hosts known because of its low redshift and intrinsic luminos- 
ity. Despite a high foreground extinction, it was already well- 
characterised in the optical, show ing a high star-for mation 
rate and relatively low metallicity. Margutti et al. ( |2007| l sug- 
gest it is somewhat atypical of the few nearby GRB hosts so 
far observed in that it lies, at low-metallicity and relatively 
high luminosity, off the luminosity-metallicity locus of KPNO 
International Spectroscopic Survey galaxies. The deviation 
might be explained simply by the low pre-existing stellar mass 
in this galaxy. Acco rding to the fu ndamenta l plane for star- 
formin g galaxies of Mannucci et al. ( 20 10); [Lara-Lopez et al.| 
( |2010[ ), HG 03 1203 should have a stellar mass of approxi- 
matel y log(M* /Mq) ~ 9. 5, similar to the value obtained by 
ICastro Ceron et dJpOlQ) . 

Initial samples of gamma-ray burst host galaxies su ggest 
they a re typically sub-l uminous and blue (|Le Floc'h et al.| 
2003 ) with young ages (Christensen et al. 20 041), high specific* 




morphologies ( ]f rucnter et al.|2006[ ). 
however have been based on rather incomplete and arbitrarily 
selected galaxy samples, affected by strong optical/UV selec- 
tion bias because of the requirement to have a detection of 
an optical afterglow to localise the burst to a high enough ac- 
curacy to identify the host. On the other hand, analyses of 
the hosts of bursts believed to be dust-obscured are consistent 
with this picture — sub-mm and radio observations yield only 
a small fraction of detections of these targeted bursts ( [Tan-| 
|vir et al.||2004[ [Berger et al. 2003) and star -formation limits 
found in the radio (Michalows ki et al.|200 8) and X-ray (Wat- 
son et al.|2004jl show no evidence of very high star-formation 



rates — suggesting that even the hosts of obscured bursts are 
not typically ultraluminous IR galaxies (ULIRGs). In partic- 
ular, the hosts of GRBs at low redshifts are consistent with 
this general outl ine (Sollerman et al. 2005 1 Michalows ki et al.] 
|2009[ |Thone et al.|20 08 ; Wiers ema et al.|2(J07) , which means 
that despite the strong evolution of star-forming galaxies from 
z > 3 to z = 0, GRBs offer a good way of targeting closer sites 
of very young star-format ion that are typ ical of galaxies in the 
early universe (e.g. Kocevski et al.|2009) . 



4.1. HG 031203 as a high-redshift star-forming galaxy 

HG 03 1203 may offer a good low-redshift archetype of 
high-redshift star-forming galaxies. Its extremely young age, 
blue optical continuum, and low stellar mass bears a striking 
resemblance to models of z and T-band drop-outs detected 
with HST. Az~2 an alogue of these galaxies was suggested 
recently by |Erb et al. (2010 1, and it is interesting to note 
that HG 03 1203 has emission line ratios and estimated stel- 
lar mass very similar to that galaxy. Furthermore, the fact 
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Figure 5. The SED of the host galaxy of GRB 031 203 . The continuous lines (blue) are the X-shooter and Spitzer-IRS spectra, open and closed squares are optical 
I Margutti et al. 2007 1 and near-IR I Prochaska et al. 2004 1 photometry respectively, while closed circles represent additional optical and near-IR photometry from 
Malesani et al. (2004) respectively. Data from Spitzer-IRAC (diamonds) and MIPS-24^tm (open diamond) as well as ATCA (open sta rs) are also plotted. Fi nally, 
upper limits are shown from MIPS-SED mode (filled, down triangles), APEX 870 /im (open, down triangles) and VLA (filled stars, Soderberg et al. 2004). The 
SED has been corrected for a Galactic extinction of E(B — V)=1.04. Marked in small open circles (grey), is the SED of the actively-star-forming blue compact 
dwarf galaxy II Zw 40 which has a similar metallicity and radiation field. Its SED is scaled by a factor of 17, the ratio of the [Ne III] 15.6 fim emission lines. The 
best-fit SED template modelled with the GRASIL code is plotted as a dotted line. 



that HG 03 1203 was selected by the presence of a GRB — 
known to be associated with the deaths of massive stars, 
and now among the most distant spectroscopically-confirmed 
objects — adds to the interest with which HG 031203 must be 
held as a nearby example of the high-redshift mode of massive 
star-formation. 

Our data on HG 03 1203 bear on the current debate sur- 
rounding the SEDs of z ~ 7 and z ~ 8 galaxies detected via 
z and T-band drop-outs. It has been suggested that the el- 
evated fluxes (above the extrapolation of the restframe UV 
continuum) redward of the 4000 A break in these objects as 
observed by IRAC in stacked 3.6 and 4.5 fim images, are in- 
dicative of a substantial 4000 A/Balmer break. From these de- 
tection s a stellar age of ~ 300 Myr is inferred for these galax- 
ies ( |Labbe et al.|2010b|at|Bouwens etaL|2010) . However it 
has been pointed out, and demonstrated via stellar synthesis 
models, that nebular line emission, in particular very strong 
[O III] lines, can prod uce an elevation in the IRAC bands sim- 
ilar to that observed (Sc haerer & de Barros|2010{ Ono et aL] 
2010} . Using models containing nebular line emission, a pop- 
ulation asyoungas5 Myr provides a reasonable fit to the data 
(Schaerer & de Barros 2010). The question of old stellar pop- 
ulations in galaxies at z ~ 7— 8 is critically important, as it has 
a strong bearing on when the first population of stars formed 
and whether these stars formed at a time that could have al- 
lowed them to cause reionisation. While it has been argued on 
a number of grounds that the elevated IRAC fluxes are indeed 
related to an old stella r population in these very early galaxies 
(Bouw ens et al.|2010| l, the arguments are circumstantial and 
the subject is still open. 

In Fig. [6] we overplot the SED from Labbe et al. ( 20 1 0b) > of 
faint z ~ i drop-out galaxies (scaled by 0.36), on the~SED of 
HG 031203. The approximate similarity in UV slope (though 
there is some uncertainty here due to the high foreground 



Galactic dust column), in the total luminosity, and in the 
relative elevation of the broadband fluxes, show clearly that 
HG 03 1203 is indeed a reasonable analogue of these z ~ 7 
candidates. Furthermore, we see in this analogue precisely 
the broadband photometry behaviour suggested by the mod- 
elling of Schaerer & de Barros (2010). The effect only occurs 
because of the stacking of sources. At z = 7, the line lies in 
neither IRAC band, but the stacking of sources in the range 
z = 6.4-7.4 and the relatively rectangular shape of the band- 
passes raises both IRAC bands. 

Lastly, in spite of the very young stellar age indicated by 
the lack of a 4000 A break, HG 03 1203 is clearly dusty, with 
as much emission from star-formation emerging in the mid- 
IR via dust-reprocessing as there is emerging in the UV. The 
best-fit GRASIL model also suggests that the UV emission is 
still significantly obscured (Ay ~ 0.3), although it is obvious 
that this must be the case from the mi d-IR emission. This 
latter p oint also echoes the claim from Schaerer & de Barros 
( |2010p based on modelling, that drop-out galaxies at z ~ 7 
may still be substantially obscured in spite of their blue con- 
tinua. From sub-mm and mm observations, we are already 
aware of the presence of large quantities of dust at z ~ 6 in 
QSOs. However, the possible presence of dust in significant 
quantities in very young galaxies at z ~ 7, suggests that dust is 
formed very rapidly. Typically, v ery rapidly is taken to mean 



a few hundr ed million years (e.g. Dwek et al.|2007 Valiante 
| et al.|[2009| ). However, if the dust is for med in association 
with the formation of star s, whether in SNe ( Dwe k et al.]2007 



Dunne et al. 2003 2009'), or grown in the ISM (Draine 2009 
Micha lowski et al.|2"0 l0c b) the timescales could in principle 
be as short as several million years. This has major impli- 
cations for the detection of galaxies of this luminosity with 
NIRCAM, NIRSPEC and MIRI onboard JWST at z ~ 7-10. 

5. CONCLUSIONS 
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Figure 6. The UV-optical SED of the host galaxy of GRB 03 1 203 (grey solid 
line and black open squares). The V-band photometry is clearly strongly ele- 
vated simply due to the [O III] emission lines. O verplotted (diamonds ) is the 
stacked SED of faint z ~ 7, z-band dropouts from Labbe et al. (2010b I scaled 
by 0.36 to match the luminosity of HG 03 1203. The blue and red diamonds 
represent redshifts of 6.4 and 7.4 respectively, indicating the approximate 
range of rest wavelengths entering the stacked filters. The elevation of pho- 
tometry near the [O III] lines and the slope of the UV-optical continuum are 
similar in both SEDs suggesting the galaxies have similar dominant stellar 
ages, < lOMyr. 

We have presented the first mid-IR spectrum of a GRB host 
galaxy as well as the first X-shooter spectrum of a GRB host. 
These spectra show strong nebular lines associated with re- 
cent star-formation, but no PAH features. We conclude that 
HG 03 1203 is an extreme star-forming galaxy with a radia- 
tion field among the hardest known for such a galaxy, in both 
its obscured and unobscured regions. The broadband SED 
of HG 03 1203 has a moderate luminosity dust emission peak, 
likely with a high temperature, similar to local blue compact 
dwarf galaxies. If dust can be formed in association with su- 
pernovae, it is possible that such a dust peak may exist even in 
very young galaxies in the early universe, making such galax- 
ies potentially much easier to detect with JWST. The lack of 
a 4000 A break and the high-ionisation emission line ratios in 
the MIR spectrum of HG 03 1203 suggests a dominant stellar 
age < lOMyr. HG 031203 has a UV-optical SED similar to 
z-band dropout galaxies at z ~ 7. This and its very young age 
indicate that this GRB host galaxy at z ~ 0.1 may be an ex- 
cellent analogue of the population of star-forming galaxies at 
very high redshifts currently being discovered by HST. 
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